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Abstract 


Consider  a  Brownian  motion  with  a  downward  drift  of  rate  a .  Its  maximum  over  all  time 
has  the  exponential  distribution  with  parameter  2 a .  Our  aim  is  to  study  this  maximum  as  a  sto¬ 
chastic  process  indexed  by  a .  That  process  is  related  to  the  convex  majorant  of  the  standard 
Brownian  motion  and,  through  the  latter,  to  a  Poisson  random  measure.  This  connection  is 
exploited  to  obtain  various  distributional  results.  The  results  are  of  interest  in  queueing  theory. 
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considered  by  NEWELL  [4]  and  by  COFFMAN,  KADOTA,  and  SHEPP  [2],  the  latter  viewing 
the  mode!  as  that  of  storage  allocation  in  computer  memory'. 

Let  Qt(n)  be  the  random  variable  that  is  0  or  1  according  as  the  stall  n  is  empty  or  occupied 
at  time  t .  The  random  vector  Q,  =  ( Q ,  (1),  Q,  (2), ...)  is  the  state  of  the  system  at  time  t .  The  pro¬ 
cess  Q,(n  );  /  £  0}  is  the  queue  size  process  in  an  M/M/°°  system;  it  is  regenerative,  and  0  is 
a  regeneration  state  for  it  It  follows  that  the  vector  (0,  0,  ...)  is  a  regeneration  state  for 
{Q,\  t  £0J  and  that  the  latter  has  an  equilibrium  distribution.  Let  Q  be  a  random  vector  (of 
zeros  and  ones)  whose  law  is  that  equilibrium  distribution. 

The  distribution  of  ^  Q  (« )  is  Poisson  with  mean  X.  The  distribution  of  £ f  Q  (n )  is  the 
equilibrium  distribution  of  the  queue  size  process  in  the  M/M/m/m  system  with  arrival  rate  X  and 
service  rate  1;  thus,  that  distribution  is  the  conditional  distribution  of  Q(n)  given  that 
X«>  Q(n)<m.  Other  than  these  facts  and  a  few  conclusions  that  can  be  drawn  from  them  by  ele¬ 
mentary'  probabilistic  considerations,  there  is  not  much  known  about  the  distribution  of  Q . 

For  a  >  0,  lei  X1'4  Y^(a,t)  be  the  number  of  empty  stalls  at  time  r  among  those  labeled 
with  n  <X-a/.-~.  ALDOUS  [1]  has  shown  that  the  process  {Y  >0:  a  >  0,  t  £  0}  converges 
weakly,  as  X  — » to  a  process  {Y  (a  ,r);  a  >  0,  t  £  0},  which  he  identified  and  showed  that,  in 
the  limit  as  t  -»  converges  weakly  to  the  process 

Y  (a )  =  map  (^2BS  -  as),  a  >  0  , 

where  B  is  the  standard  Brownian  motion.  He  calls  fY (a);  a  >0}  the  exponential  process,  after 
the  well-known  fact  that  Y{a )  has  the  exponential  distribution  with  mean  1/a  for  each  a . 

Our  main  contribution  is  to  supply  the  probability  law  of  Y  in  simpler  terms.  For  this  pur¬ 
pose  we  choose  to  work  with 

(1.1)  Za  =  -LrC'fJa)  =  map  (B,- at),  a>  0, 

and  let  Da  be  the  last  time  r  at  which  B,  -  at  touches  its  zenith  Za ,  that  is, 
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(1.2) 


Da  =  sup  {t  :  B,  -  at  =  Za)  ,  a  >  0  . 


It  turns  out  that  Da  is  the  left-derivative  of  Z  at  a  and,  thus,  is  related  to  the  density  of  empty 
stalls  in  the  parking  lot,  in  equilibrium,  around  k  -  a  k3/<  for  large  k. 

The  next  section  contains  a  few  simple  geometric  observations.  First  we  relate  the  process 
(D ,  Z)  to  the  convex  majorant  of  the  Brownian  motion  B  .  Using  the  hard  results  of  GROENE- 
BOOM  [3]  and  PITMAN  [5]  about  the  latter,  we  are  able  to  express  D  and  Z  in  terms  of  a  Pois¬ 
son  random  measure  on  (0,»)  x  (0,<»).  It  follows,  in  particular,  that  D  has  non-stationary 
independent  increments,  and  that  (D  Z)  is  a  non-homogeneous  Markov  process. 

The  process  a  — >  Za  is  continuous,  concave,  and  decreases  from  its  limit  +°°  at  a  =  0+  to  its 
limit  0  at  Therefore,  its  "hitting  time"  process 


(1.3) 


Az  =  inf  {a  :  Za  <z}  ,  z  >  0 , 


is  the  functional  inverse  of  Z .  It  turns  out  that  the  process  A  has  the  same  probability  law  as  Z . 
This  observation  is  also  put  in  the  next  section. 

The  last  section  is  devoted  to  computational  issues.  We  compute  the  joint  distribution  of 
Da  and  Za  and  also  the  transition  function  of  the  Maikov  process  (£>  X). 


2.  ZENITH  PROCESS 


The  problem  with  the  definition  (1.1)  of  Za  is  that  it  suggests  re-drawing  the  path 
t  — »  B,  -  at  if  we  wish  to  van1  a .  The  following  observation  circumvents  the  problem: 


(2.1)  Za  =  inf  {x  >  0  :  x  +  at  >  B,  for  all  t  k  0}  . 


Obviously,  this  is  a  re-wording  of  (1.1),  but  the  mental  picture  it  suggests  is  much  more  con¬ 
venient  for  manipulating  a :  the  line  t  -*  Za+ at  is  the  infimum  of  all  lines  of  slope  a  that  never 
touch  B .  This  picture  is  drawn  in  Figure  1  below. 
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Let  C  denote  the  convex  majorant  of  B  ,  that  is,  the  minimal  convex  path  that  dominates  B 
(sec  Figure  1  again).  Note  that  Za  and  Da  are  determined  by  C :  the  line  t  -» Za  +  at  is  the 
infimum  of  all  lines  of  slope  a  that  never  touch  C ,  and  Da  is  the  last  time  at  which  that  line 
touches  C .  In  fact,  for  fixed  a  >  0,  Da  is  almost  surely  the  only  time  t  with  C,  =  Za  +  at . 

It  is  known  (see  GROENEBOOM  [3]  for  instance)  that  C  is  continuous  and  piecewise 
linear.  The  countable  collection  of  its  vertices  has,  almost  surely,  only  one  accumulation  point, 
namely  (0,0).  Fix  an  a  >  0;  note  that  (D„,  Za  +  aDa)  is  a  vertex;  let  To,  T  i,...  be  the  lengths  of 
successive  intervals  of  linearity  going  to  the  right  from  Da ;  let  7_i,  7"_2,...  be  those  to  the  left; 
and  let  S,  be  the  slope  of  C  over  the  interval  whose  length  is  denoted  by  7).  The  following  major 
result  was  obtained  by  GROENEBOOM  [3];  a  simpler  proof  using  the  excursions  of  B  may  be 
found  in  PITMAN  [5J. 


(2.2)  THEOREM.  The  pairs  (S, ,  7) ),  i  e  Z,  form  a  Poisson  random  measure  N  on 
(0,°°)  x  (0,oo)  whose  mean  measure  is 


(2.3)  v(ds,dt)  =  £.  yM), 

where  ys  is  the  gamma  distribution  with  shape  index  1/2  and  scale  parameter  s2/ 2  (the 
corresponding  mean  is  1/s2). 

The  probability  law  of  a  Poisson  random  measure  is  determined  by  its  mean  measure. 
Thus,  the  following  specifies  the  probability  law  of  ( D,Z ).  For  computational  purposes,  the 
representations  given  here  for  Da  and  Za  are  the  key  starting  points. 


(2.4)  PROPOSITION.  For  each  a  >  0, 


Da  =  J  N(ds,dt)t  , 

la  ,«o}x(0,«») 

oo 

Za  =  f  Dsds  =  f  N(ds  ,dt)  (s  -a)t 

a  [a  ,«»Jx(0,oo) 


K* 


A 
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The  process  D  has  non-stationary  independent  increments.  The  process  (D ,  2 )  is  a  temporal]) 
non-homogeneous  Markov  process. 


PROOF.  First  note  that  (see  Figure  1) 


Da  =  X  Ti  1  (a  ,oo)  (S,) , 


Za  +  a  Da  —  B  (Da)  —  £  Si  Ti  1  [a ,«>)  (Si) 


Expressed  in  terms  of  the  Poisson  random  measure  N ,  these  become  (2.5)  and  (2.6).  The  remain¬ 
ing  statements  are  immediate  from  the  independence  of  the  restrictions  of  N  to  disjoint  Borel 


Figure  2  below  shows  the  qualitative  features  of  D:  it  is  piecewise  constant,  left- 
continuous,  and  decreases  from  its  limit  +°°  at  a  =  Of  to  its  limit  0  at  +°°.  It  follows  from  (2.6) 
that  Z  is  continuous,  concave,  piecewise  linear,  and  decreases  from  its  limit  +°°  at  a  =  0+  to  0  at 


FIGURE  2 
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It  was  noted  by  ALDOUS  [1]  that,  for  each  c  >  0,  the  process  cZca)a>c  has  the 

same  probability  law  as  ( D,Z ).  This  can  be  seen  from  the  preceding  characterization: 
a  — »  c2Dca  jumps  at  the  points  St/c  by  the  amounts  c2Ti\  the  pairs  (S./c,  c2Ti)  form  a  Poisson 
random  measure  that  has  the  same  mean  measure  as  N\  hence,  a  — >  c2Dca  has  the  same  law  as 
D. 

We  end  this  section  with  an  observation  on  the  process 

(2.7)  A2  =  inf  {a  :  Za  <  z }  ,  z  >  0  . 

Obviously,  z  — >  A:  is  the  functional  inverse  of  the  one-to-one  mapping  a  —*Za  of  (0 ,«>)  onto 
(0,oo).  u  follows  that  the  qualitative  picture  of  A  is  exactly  that  of  Z.  In  particular,  A  is  piece- 
wise  linear  and 


(2.8; 


D2  =  lim 
c*o 


Aj  +c  —  A: 


1 


z  >0. 


e  ~&lAj  * 

The  process  D  is  piecewise  constant,  right-continuous,  decreasing. 

(2.9)  PROPOSITION.  The  process  (D ,A)  has  the  same  probability  law  as  the  process  (D , Z) 
In  particular,  the  collection  { Z(S,) ;  i  €  Z]  has  the  same  law  as  the  collection  {S,\  i  eZj;  they 
form  Poisson  random  measures  on  (0,«)  with  mean  measure  dsls . 


n; 

ft 


j' 

a 


PROOF.  We  put  (2.7)  and  (2.1)  together  and  manipulate: 

Az  =  inf  { a  :  inf  {x  :  x  +  at  >  B,  for  all  t)  <  z} 
=  inf  {a  :  z  +  at  >  B,  for  all  t) 


\ 


the  zenith  process  associated  with  B  .  Since  ( uB  1/K)  is  a  standard  Brownian  motion  like'  B ,  it  fol¬ 
lows  that  A  has  the  same  probability  law  as  Z .  This  proves  the  first  statement,  since  D  is  the 
derivative  of  -A  and  D  is  the  derivative  of  -2 . 


The  points  S,  are  the  jump  locations  of  D ,  and  the  points  Z  (S, )  are  those  of  D .  This  proves 
the  second  statement. 


r 


3.  ENTRANCE  LAW  AND  TRANSITION  FUNCTION 


We  derive  the  distribution  of  the  random  variable  ( Da ,  Za)  and  the  transition  function  of  the 
process  (D  ,Z).  The  computations  rest  on  the  characterization  given  by  Proposition  (2.4)  and  on 
the  well-known  formula  for  the  Laplace  functional  of  the  Poisson  random  measure  N  with  mean 
measure  v: 


(3.1)  £  exp  -|N(dlr)/(x)  =  exp  -  J  v(dx )  (1  -  e  -OO) 
for  every  positive  Borel  function  /  on  (0,«=)  x  (0,*>). 

(3.2)  PROPOSITION.  For  each  a  >0, 


(3.3)  £  exp  (,-pDa -qZa)  =  - - - ,  p>0,  q>0; 

a  +  q  +  ^a2  +  2p 

(3.4)  P{Da  e  dt,  Zae  dz  }=dtdz  2az  eJ^= f  r  >  0,  z  >  0  . 

V27tf3 

In  particular, 

#6  (J 

(3.5)  P{Za  €  dz}  =dz 2ae ~2ni  ,  P  {Da  e  dt}  =dt  \  du  ~=  . 

1  <2nu3 

PROOF.  Fix  a  >  0,  p  >  0,  q  £  0.  In  view  of  (2.5)  and  (2.6), 


pDa  +  qZa  =  J  N(ds  ,dt)  (pt  +  q (s  -  a)t) . 

Using  (3.1)  and  the  form  of  the  mean  measure  v  given  by  (2.3),  the  Laplace  transform  (3.3)  is 
obtained  via  elementary  calculus.  To  invert  the  Laplace  transform,  first  write  it  as 
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|  dz  e -*z  2 ae~*:  e-*^ai  +  2P 

and  then  recall  that  e~2^ ,  r  20,  is  the  Laplace  transform  of  Hz,  the  first  time  a  standard 
Brownian  motion  hits  the  level  z ,  that  is, 

e-2<a>  +  lp  _  7  dt  .  z  e-ja/Zf  eApMxVl)i 

The  rest  is  trivial. 

(3.6)  REMARK.  Although  (3.4)  is  explicit  and  shades  of  exponential  and  stable  distributions 
can  be  felt,  it  does  not  seem  well-suited  for  probabilistic  thinking.  The  following  representation 
is  better,  especially  for  Monte-Carlo  methods.  For  a  >  0, 

a2Da  =  X  (1  ~<U)2,  dZa  =  X<U  (1  -<U), 

where  X  and  U  are  independent,  U  has  the  uniform  distribution  on  (0,1),  and  X  has  the  gamma 
distribution  with  shape  index  3/2  and  scale  parameter  1/2. 

The  following  specifies  the  joint  Laplace  transform  of  any  number  of  increments  of  Z  (upon 
taking  f  =  p\  1A]  -K..+  p„  1A>  with  A  An  disjoint  intervals). 

(3.7)  PROPOSITION.  For  any  positive  Borel  function  /  on  (0,°o), 

oo 

E  exp  J  /  (a )  dZa  =  exp  -  f  ds  (1  -  .  ■  ■  _  ) 

W-)  t)  s  Vs2  +  2/(s) 

where  /  (s )  is  the  Lebesgue  integral  of  /  over  (0  j ). 


PROOF.  Note  that 


\f(a)dZa  =  -j  f(a)Dada  =-j  N(ds,di)  f(s):  , 

and  use  (3.1). 

As  mentioned  in  Proposition  (2.4),  the  process  D  has  non -stationary  independent  incre¬ 
ments,  and  the  process  (D  Z )  is  a  non-homogeneous  parameter  Markov  process.  Let 

(3.8)  Pab(tJ  ;  du ,  dy )  =  P{Dbe  du,Zbe  dy  I  Da  =t,Za  =  x) 

for  0  <  b  <  a ,  0  <  t  <  u ,  0  <  x  <  >  (in  our  zeal  to  deal  with  positive  random  variables,  we 
choose  to  work  with  the  parameters  in  decreasing  order).  In  view  of  (2.5)  and  (2.6), 

(3.9)  Pab(tjc  ;  du,dy)  =  P (t  +U  e  du  ,  x  +(a  -b)t  +Y  e  dy)  , 
where 

(3.10)  U  =  f  N(ds.dt)t,  Y=  [  N  (ds ,  dt )  ( s-b)t 

(b  fi  N(0,«=)  |b^)x(0,">) 

The  joint  Laplace  transform  of  U  and  Y  can  be  obtained  from  (3.1)  as  in  the  first  step  of  the 
proof  of  (3.2): 

(3  11)  E  e~Pv-*Y  =  —  •  -  +C}  +  2  4  4  2(a  ~b)q 

a  b+q  +<b2  +  2p 

=  b_  +  a  -fr  .  1  .  2b 

a  a  1  b+q  +V^2  +  2p 

+  id  '  - 2h _ i  J  fr4q)<fr 

c  2  b  +  <?  +'lb2  +  2p  a~b  *  V{c  +q)2  +  2p  -2bq 

Inverting  this  is  tedious  but  manageable.  It  yields  the  following  for  the  distribution  of  the  pair 


i*4w»**. 


(U,Y): 


r*rr 

E? 


(3.12)  cp  =  L  6(0  0)  +  (1  -  £)  kb  +  \  kb*  ) 


where  the  asterisque  denotes  convolution,  hx  is  the  Dirac  measure  at  x ,  kb  is  the  distribution  of 
( DbZb )  specified  by  (3.4),  and 


(3.13) 


\iab(dt  ,dy)~  -d'  h(a-t,)t  (<#>■).  r  >  0,  y  >  0  . 

V2ru3 


Putting  the  distribution  of  (U  ,Y)  into  (3.9)  yields  an  explicit  expression  for  the  transition 
function  Pab .  As  a  by-product,  we  have  the  joint  distribution  of 


U  =  Db  -  Da,  Y  =  Zb  -  Z0  -  (b  -a)Da  . 

Noting  that  Da  is  independent  of  {U  ,Y),  one  can  obtain  the  distribution  of  (Db  -  Da  ,  Zb  -Za  ) 
among  other  things. 

However,  it  is  clear  that  such  results  are  of  limited  use  because  of  their  complexity. 
Overall,  the  computational  complexity  is  caused  by  a  confluence  of  two  incompatible  operations, 
addition  and  multiplication:  look  at  the  form  (2.3)  of  the  mean  measure  v;  the  Haar  measure  ds  Is 
indicates  that  the  natural  group  operation  on  the  jump  points  S,  is  multiplication,  whereas  the 
jump  amounts  T,  are  additive. 

Of  course,  it  is  easy  to  transform  the  Poisson  random  measure  N  into  one  with  a  nicer  inten¬ 
sity:  define  /  to  be  the  mapping  (j  ,t )  ->  ( log  s,  s2t );  then  the  image  of  A'  under  /  is  the  Pois¬ 
son  random  measure  N  =  A '/  _1  on  (-»,«>)  x  (0,«>)  with  mean  measure  duy(dv)  where  y  is  the 
gamma  distribution  with  shape  and  scale  parameters  equal  to  1/2.  But,  then,  expressions  for  D„ 
and  Za  in  terms  of  A  have  to  undo  the  transformation,  and  there  is  no  gain  at  the  end.  Using 
p  -  log  a  to  index  the  processes  involved  (and  working  with  Zp  =Z(ea))  does  not  help  either. 


On  the  other  hand,  it  is  easy  to  describe  the  construction  of  the  path  of  (D  Z)  over  the  inter¬ 
val  (0 ,a  ].  This  may  be  useful  for  Monte-Carlo  purposes. 

First,  we  observe  that  the  conditional  distribution  of  S,+i  given  (Si,  S,_i,...)  is  the  uniform 
distribution  on  (OS,)-  Thus,  to  construct  (D  Z)  over  (0 ,a  ],  we  start  with  U  and  X  described  in 
Remark  (3.6)  and  generate  D0  and  Za-  Then,  we  let  U  i,  U  2,...  be  i.i.d.  uniform  on  (0,1),  let 
X 1,  X  2»...  be  i.i.d.  Gaussian  with  mean  0  and  variance  1,  and  set 

(3.14)  Sc  =  a  ,  S,  =  aU  \  U 2  ' ' '  V i  ,  T,=(-^-)2,  i  =  1.2 . 

With  these,  define 

(3. 15)  D(Sa)  =  Da  ,  D  (S,)  =  D  (S^)  +  T,  .  i  >  1  , 

(3.16)  Z(S0)=Za  ,  Z(S,)  =  Z(S,_i)  +  (S,_i  -  Si)  •  D(S,_i) ,  i*l. 

Then,  ( Db ^t(0<1  is  the  left-commuous  piecewise  constant  path  whose  value  at  S,  is  D  (S, ),  and 
(Zfe).«(o,j  *s  the  continuous  piecewise  linear  path  whose  value  at  S,  is  Z(St).  Incidentally, 
(S, ),  (Si ,  D  (Si )) ,  (Si ,D(S,),Z (S, ))  are  all  Markov  chains. 
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